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Abstract 

The stationary condition (Hopf equation) for the (n+1) point correlation function 
of a passive scalar advected by turbulent flow is argued to have an approximate SL{n, R) 
symmetry which provides a starting point for the perturbative treatment of less symmetric 
terms. The large scale anisotropy is found to be a relevant fleld, in contradiction with 
Kolmogorov phenomenology, but in agreement with the large scalar skewness observed in 
shear flows. Exponents are not universal, yet quantitative predictions for experiments to 
test the SL{n, R) symmetry can be formulated in terms of the correlation functions. 
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Kolmogorov succeeded, by very simple arguments (K41), in predicting the scaling 
for the velocity correlations < VkV-k > at high Reynolds number R, for wavenumbers 
intermediate between those defined by the geometry of the flow and dissipation; the so 
called inertial range. ^ Obukhov and Corrsin soon applied the K41 reasoning to a passive 
scalar i.e., a field obeying, 

dtO + v-Ve = kV^Q (1) 

(k is diffusivity), and experiments subsequently found scaling behavior for Q, although 
the K41 exponent is approached only at very high R if at all.^ However an even more 
glaring inconsistency with K41 appeared with the observation that the derivative skewness 
Sd =< {djf > / < {djf >^ is of order one and R independent out to highest R 
available.^ This fact motivated our work. 

Since Sd breaks parity we follow the experiments and impose a large scale gradient g so 
that Q = 9+g-r where < 9 >= and < 6'(r)6'(0) > inherits a correlation length or integral 
scale from the velocity field. The inhomgenous term g ■ v which appears in the $ equation 
then acts as a "force" which maintains the 9 fluctuations stationary. The conundrum with 
Sd is that K41 predicts small scale isotropy for large R specifically, Sd ~ 

or for 69r = 9{r) - 9{0), S =< 69^ >~ ri, (vs in experiments). The K41 
arguments have proved so seductive and work so well for the 2-point velocity correlations 
that problems with the scalar and difficulties with higher order velocity derivative statistics 
in simulations have been all but forgotten. 

The observation of Sd ~ 0{1) is particularly intriguing since it suggests that ^ is a 
relevant variable along with the "energy" (scalar variance really) dissipation rate e, in 
the sense of parametrizing the effect of the large, geometry specific scales on the inertial 
range. Averaging g away by insisting on isotropic large scales will eliminate s^, but in no 
way lessen the instantaneous non local interactions between large and inertial range scales. 
Whether this merely obscures the relevant physics, much like studying superfluid Helium 
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in a fixed number ensemble, or actually eliminates the dominant exponent for even order 
correlation functions remains to be seen. 

The prospects for an analytic theory of scalar intermittency are much enhanced by the 
observation that remains ~ 1 and R independent even when the complex turbulent shear 
flow is replaced by a Gaussian v fleld.^ Other simulations and a plausible closure scheme 
of Kraichnan make the same point for even powers of S9r^: boring but scale invariant 
velocity fields generate interesting scalar statistics. 

In this report we argue that the dominant term in the evolution equation (the "Hopf 
equation" ) of the scalar multipoint correlator i/j generated by a velocity field with a phys- 
ical correlation time is highly symmetric and is integrable by Lie Algebraic methods. The 
leading (anomalous scaling) terms in ifj arise as the zero modes of the Hopf operator, 
Lh, with the lowest scaling index A.^'^'^ For the dominant, symmetric part of Lh, the 
zero modes are labeled by a complete set of quantum numbers, but because of the high 
symmetry the A are infinitely degenerate, i.e., independent of a subset {q} of the quan- 
tum numbers. Smaller, lower symmetry, perturbations, which are nevertheless dominant 
when points are coUinear are treated by singular perturbation theory, or equivalently by 
diagonalizing within the degenerate subspace. Like the operator itself, the exponents and 
correlator depend on the details of the velocity field beyond its scaling dimension. However 
since with a general configuration of points is governed approximately by the symmetric 
part of Lff, it can be represented as a linear superposition of the degenerate {q} modes 
with the weight defined by i/j with points coUinear. The latter can be measured directly in 
experiment or numerical simulations and the above relation used as a quantitative check 
for the existence of the SL{n, R) symmetry. 

The homogeneous part of the inertial range Hopf equation for a white noise velocity 
field scaling as < 5v^ >~ r'^ may be written in terms of the Richardson operator 

L^R^ = r-i"^ {s-\d + 1 - C) - (2 - CKrtj)dtd^ (2) 
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acting on all pairs of points {i,j,fij = Vi — fj) and where a,b denote spacial indices.^ 
However within K41 theory (C = |) the velocity is not white but has a scale dependent 

2 — 1 

(Lagrangian) correlation time ~ r'Se~ . This is just a restatement of the scale invariance 
of the velocity equation which says that the change in spacing r of two points over the 
correlation time of their relative velocity, (5fr, is ~ r. Take the spacial arguments of some 
multipoint correlation roughly equidistant, let Rg be of order their radius of gyration, then 
a first approximation to the action of the velocity field is an 0(1) volume preserving linear 
mapping of all points. For Gaussian v, the averaged operator implementing this finite 

(2) 

coordinate change is just the exponential of Lq = Lj^ , so the condition of stationarity can 
be expressed as Loip = 0. (Note Lq is identical with the Batchelor-Kraichnan operator^° 
for white random gradient advection but its meaning here is quite different.) The action 
of V from scales less than Rg is incoherent and uncorrelated between points and plausibly 
white in comparison with the coherent part. Upon averaging it gives rise to an "eddy 
dissipation" operator Ld (acting along with Lq) the form of which we will discuss later. 
Finally, the velocity on scales ^ Rg is an overall translation which in the presence of the 
large scale gradient generates an inhomogeneous forcing term, which however does not 
affect the zero modes. 

The operator Lq is a very attractive starting point for perturbation theory since 
it is integrable. For the (n + 1) point correlation function of eliminate the center 
of mass by defining n difference vectors, pi — (ri — r^)!^^-, P2 = (^'i + r2 — 2r3)/\/6, 
P3 = (^'i + r2 + rs - 3r4)/ \/l2, etc. 

Lo(n) = -{d + 1)L^ + 2dG'^ + d{d - n) (A^/{nd) + a) (3) 

where L^ = ^^{pfd^ — p\dfY is the square of the total angular momentum; A = p^df, 
the dilatation operator; and = ^^GijGji, Gij = p^dj — ^5ijA, is the Casimir of 
the group SL{n, R) acting on the i— label, diagonal in space, henceforth called isospin. 
The appearance of the SL{n, R) in (3) is the consequence of the underlying evolution step 
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being the multiplication of all pi by a spatial strain matrix — the linear mapping described 
earlier — which is invariant with respect to a general linear transformation acting on p's. 
One might have expected the group SL{d, R) (Casimir J^) to appear, since we are applying 
unimodular matrices in space, and it does if we exploit a "duality" relation on to replace 
the last 3 terms in Lq by 2dJ'^. 

The zero modes of Lq are found by constructing the representations of SL{n)xSO{d) x 
A. The SL groups are non compact, their representations infinite dimensional, and those 
of interest non unitary — we do not have a suitable inner product. We work in the space of 
homogeneous functions whose scaling dimension. A, diagonalizes A. The remaining quan- 
tum numbers are discrete. As before we define v, related to A/2 by an integer, (defined 
below) in terms of which ^ (n — l){2iy/n + (for n < d). These quantum numbers, 
together with those for SO{n) (which is one choice for the maximal commuting sub algebra 
of SL{n)) and SO{d) are a maximal set for (n, d) = (2,2) and the corresponding eigen- 
functions were given in [9]. The most relevant (i.e., smallest non negative[ll]) eigenvalues 
A for any (n, d) are A = 1 for p - wave, i = 1; and A = for s - wave. The latter could 
have been anticipated from the Batchelor-Kraichnan interpretation of (3). Finally we note 
that acting on any function behaving as a; > the operator Lo(?i) reduces in the 

limit Tij — > smoothly to Lo(n — 1) acting on the remaining variables. This is because 
Lq originates from the advective part of (1) which evolves d'^{r) the same way as ^(r), 
allowing the different points in the correlator to be brought together with impunity. 

The dissipation in (1), and its restatement in terms of Ld, the incoherent eddy damp- 
ing part of the Hopf operator, of course can not reduce so simply. It does matter when 

(1) - (-) 

two points are brought together. The most plausible form of is L^-^^ = aRgL^ via 
(2) where i?| = (pf)3 accounts for the time scale implicit in Lq. Alternatively, i?| just 
gives to Lo the same scaling dimension as Lq; there are no intrinsic scales in the inertial 
range. We will do perturbation theory in the parameter a, the structure of which convinces 
us that it is equally plausible to use an analytically more tractable model which induces 
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the same local singularities into the eigenfunctions ip, L\y = aRg P{p)'V^ where /?(p) is 
homogeneous of degree |, inherits the permutation symmetry induced by rj on the 
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p's, and P{p) ~ r?j as any pair coalesce; a "Kolmogorov point". Near such a point 
ip has local s-wave symmetry in r^. It will also be useful to consider a scale invariant 
Laplacian model = (^)-^a^p- 

For CK <S 1, L£) is a singular perturbation to Lq, in particular near a Kolmogorov 
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point we can chose the pi such that pi — > and Ld ^ Pi^li whereas Lq involves powers 
of Pidp^. Now when pi — > the homogeneous Hopf equation becomes, 

4^"'^ < ^2(1)^(2).. > +aRl < e{l)e{2).. >= (4) 

(plus 0{a) less singular terms in the case of L^^^). The new operator e is the local dis- 

(-) 

sipation rate, expressible as L)^ {p)9{p + 1)^(1) in the inertial range and matching to 
k(V6')^ in the dissipation range. Of course (4) is not a closed equation, but < e9.. > 
must have a positive scaling dimension, (which can be demonstrated by observing from 
(1) that sup|^(r)| and < e > are finite and set by the large scales). Hence, the exponent 

r 

for < 66... > must exceed |.^^ Scaling and the generic behavior of the correlator when any 
two points coalesce, insure that e carries a dimension of | relative to 6'^ in all correlators. 

To proceed with an explicit calculation, the physical coordinates p can be factored, e.g. 
in the case (n, d) = (2, 3): p" = 'Rmix) Vi/ where R represents isospin rotations by x 
and fji^2 are two orthogonal unit vectors defining a representation for SO (3); {fj3 = f}\/\f}2)- 
Suitable, "gauge" conditions must be imposed on functions defined on the Xi C) V manifold 
if they are to reassemble into smooth functions of p. 

Let h^{xoi) be a homogeneous function of degree v in its arguments, a = 1 — 3, then 
any function of the form ip^ = j e^'^'^h^{nip1e'^)D\^^^{e)d(l)de is an n = 2 eigenvector of 
and with eigenvalues v{y -\- 1) and 1(1 -\- d — 2) provided ni + in2 — e**^ and D is a 
representation of /SO (3) defined by the triad Cq,. We can also multiply ij^o by any element 
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of Pi A p2 since GijPi A p2 = 0. The p dependence of all these functions can be written 
interms of (w = 2Ci + Cl)) 



which has a complete set of 6 quantum numbers and defines the sense in which (3) is 
integrable. An integral representation for the Jacobi function, P, was already given in Eq. 
(6) of Ref. (6), and v is related to ^ by an integer due to boundary conditions at w = 0. 
The index m' corresponds to rotating on the isospin index of f}i^2 and is not summed. 
Acting on a general function of r]i) with the A, and Xi dependence of as in (5) 

where 1^ = i~^{r]i ■ 82 — r]2 ■ di). Reexpressing the Laplacian, 



e-^''^(66)^ -i^ ViV' = k(l - w^)d^ + (1 + A)(l - ^2)^-1^^ + -rXVw^+ 

For the skewness, the most relevant eigenvalue of (3) has Z = 1 so we can write 
(f) = (/)i{w)r]i + i(/)2i'w)r]2, then a^; (cr are the Pauli matrices) acting on ((/>!, (/)2), 

(773 • da)^ r]i, = —Sab, and |^ is a function of w. It suffices to treat < w < 1, negative w 
are equivalent to a similarity transformation of the equation by u^. 

Thus the Hopf operator with Laplacian damping, Lq-{-Uj^, reduces to a pair of second 
order equations, for which the eigenvalue problem can be solved numerically by shooting 
from near ~ 1, where Lq dominates; to ~ or pi^2 parallel {n,b,pi A p2 = CiCa^a) 
which is an invariant set under the coherent part of the velocity field. The eddy damping, 
i.e., (6b), dominates there and our normalizations require us to impose that (f) diverge 
no more rapidly than w~ . The most relevant eigenvalue with the 3-fold permutation 
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symmetry required for < 699 > is q = ±3, followed by g = ±9, since there is a gauge 
requirement for odd q when / is odd. We find A — 1 = — 3.066q;5 + 1.436q!^, the 
first term is analytic (see below) while the others are fit to numerical data accurate for 
a <1. Thus A is nonuniversal, i.e., a dependent, but stays close to its anomalous value in 
experiments of ~ 1, all the more so since the correction term need not be positive (A3 need 
not be greater than 1); e.g., for certain models of K41 eddy damping. The eigenvector can 
be thought of a superposition of modes of the form (5) with respect to the discrete allowed 
V at fixed A.^ 

We now set up perturbation theory via matched asymptotic expansions with sufficient 
generality that it applies to the L)-)' case with j3 = Rg j3 /4, j3= Po{x)Po{x~ f )/^o(x+ f )> 
(3o = 2-3(1 - cos(2x)Vl - ^2)3. The Kolmo gorov points correspond to w = and 
X = 0, and we have made the convention < ^2- The transformation x ^ x ± ^ 

corresponds to a permutation and reflection of rj implying that is antiperiodic in x with 
period ^ and of course periodic with period Also (f) is even under interchange of ri,2 
i.e., pi -pi p2 P2 or (x, Ci) ^ -(x, Ci), 6)^1 unchanged. For w = 0, is even 
around x = 0) which together with its antiperiodicity implies 0(x = ±-|) = O. 

To go through the region ~ a where Lq and Ld are comparable, we have to solve 
the following crossover equation, (0= w~^(f)) 



-^(l + ^.)/«'')5+^(A-l) 0=0 

where we have replaced q^ by —d^. Near the Kolmogorov points generally (3^ (w^ + 4x^)^ 
and the 0{a) derivative terms in (7) will be recognized as a 3-D Laplacian in cylindrical 
coordinates {w,x)- Elsewhere <^ dw For a <^ w'^ < 1 the solution of (7) must reduce 
to a linear combination of solutions (5) to (3) with u = 1/2 and scaled with Rg which for 
small w behave as 0° =^||g|t(; + |(1 — q^)w'^, (1 — ^q^w^)sgn.{q)j / {q'^ — 1). This has to 
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be propagated via (7) to < <^ a which is particularly simple for 02 which is found 
to behave as (a(x) + ^d^a w'^ ...). Substituting into (7) we derive the eigenvalue problem 
along = as 

~ 20 
a/?(x)(aJ + l)a + -(A-l)a = (8) 

The eigenvalue for Laplacian dissipation follows by setting (3= 1 and a = e^'^^ while in the 
Kolmogorov case we integrate from x = to the antiperiodic point X = f with boundary 
conditions a ~ (^1 + 0(1 - \)\x\^ + 0(x^)) and find (A -!)/«= -0.0707 (versus -0.2 
for d=2). The eigenfunction is monotonic and largely fixed by the symmetries and the 
boundary conditions at x = 0, ^. 

Neither the computed a(x) nor A can be directly compared with the experiment since 
they depend on details of the eddy damping. However we can make testable predictions 
by expressing =< ^i^2^3 > for Pi parallel (^i = by convention and tan(x) = P1/P2) 
as ^2^(x)^2 • g to within an overall constant. Then away from colinearity {uP' > a), the 
solution is a superposition of "SL(n) modes" (5) weighted by a^, the Fourier transform of 
a(x)- More explicitely we use the integral representation of and write a convolution 

Jo (cos2 ((/.)e? + sm^ i(t>)$l)^-^ 
where hq = sgn(g)(g^ — l)aq defines the q values in (5) that are superimposed, u — ^, 
and we summed on mf to impose the symmetry under ri,2 interchange. To be consistent, 
we should expand the denominator to first order in (A — 1) (the 0^*^ order term vanishes). 
Then for w <C 1, 0° is just the transform of the (^,?7) factor in (9) which explains the 
relation of h to a. 

The fiatness in two spacial dimensions is illustrative of correlation functions with 
n > d. By interchanging the meaning of space and isospin labels we can use the same X) 
^, rj coordinates as before but to emphasize their new meaning we put quantum numbers 
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p, q, qi on D in (5) and replace by m — > for the s-wave state. Equations (6a, h) also 
remain valid though it is now convenient to recognize that the 773 ■ da (a=l,2) together 
with /a comprise a spin algebra operating on the qi index of Dqq,{r]). 

For n = 2, the isospin was Abelian, and for Laplacian damping the most relevant 

eigenmode had g = 3. The quantum numbers for the analogous flatness mode follow by 

4 

writing JJ^^r • interms of pi which gives a particular combination of {q,p) = (0, 4) and (3, 

1 

4). Note only p and qf appear in (6b) not q. Thus the solution ~ '^^cp'^' {w)Dqg,{r]) to the 

q' 

Hopf equation Lq + = can be sought in a 2p + 1 = 9 dimensional space labeled by 
qf which further symmetries reduce to a three dimensional sub space (0, 2, 4) even under 
q/ ^ -qf. The eigenvalue is Xp = %y/2a/S + 2.2a - 3.8a^/^ + 1.64a^.... As before the 
leading term can be obtained analytically by solving the crossover equation. 

The analytic details of the perturbative theory for A f and its eigenmode for either L^-* 
or damping resemble the n = 2 case in that a solution to (3) ioi a <^ < 1 has to 
be matched to one for ^ <^ a. The geometry differs because before the Kolmogorov 
points, e.g. |pi| = or x = = 0, were subsets of the condition for parallism, w = 0, 
whereas for the 2D flatness the Kolmogorov points are for instance r]f = 0, (a = 1, 2 is 
space), while the condition all three pi parallel is still w = 0. The two sets have only 
a lower dimensional intersection. The eigenvalue problem analogous to (8), still operates 
within w = but the space is S0{3). The integral representation replacing (9) will apply 
for ^ a and also away from the Kolmogorov points. 

It is noteworthy that the dissipation is a bigger perturbation on the flatness than the 
skewness i.e., y/a vs a; but we have not seen how A_f > | emerges from the calculation, 
so for the moment we take this as a lower bound on a. 

The symmetry breaking external gradient ^ is a relevant pertubation in that normal- 
ized high order odd moments of d9, which vanish under isotropic conditions, will grow 
as R^' when pn = N/3 — Xn < where A at is the N point d correlator. Thus the 
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derivative skewness may already diverge (A5 < 1 for the pseudo Kolmogorov dissipation 

), and for small a all higher odd moments certainly diverge, since Ajv = 1 + 0{a)}'^ 
Note for even N, ji^ is also the scaling exponent for the N/2 point dissipation correlator, 
i.e., < e(l)e(2).. >~ R-t"^ . 

Finally, a symmetry based classification also leads one to ask qualitative questions such 
as whether for physical values of the dissipation the nominally subleading p-wave exponent 
for flatness and other even moments could be more relevant than the s-wave exponent and 
whether there could be a scaling crossover as the points depart from coUinearity. 

Our phenomenological approach is based on an approximate symmetry and its experi- 
mental consequences, e.g. (9), should be more general than the specific equation we solved. 
We hope also that the SL invariance is a more robust property of the dynamics than the 
scaling exponents and therefore visible at lower Reynolds' numbers. Our dissipation was 
adhoc, and for n > 3 should look more hierarchical. Yet compared to the study of the 
exact Hopf equation for the white noise^'^'~^° exploring different dissipation models has 
the advantage of confronting one with the non-uniqueness and the non-universality of the 
problem at hand. In particular, exponents are non-universal.-^^ 

This conclusion may be regarded as a peculiarity of the passive scalar problem defined 
here for arbitrary random velocity fields, in contrast to the properly turbulent velocity 
field governed by the Navier-Stokes equation which for large R, according to K41, posseses 
universal statistical properties. Yet, our analysis illustrates how the large scale anisotropy 
can affect the small scale statistics. Applied by analogy to the velocity field, it would 
suggest the possibility that the small scale velocity statistics may in fact depend on the 
large scale boundary conditions leading to non-universal exponents for sufl&ciently high 
order correlators.^^ This of course would just be a manifestation of the perhaps not so rare 
"rare" events directly coupling the large and the small scale of the fiow. 
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